p-meson, Bethe-Salpeter equation, and the far infrared 
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The Bethe-Salpeter equation in QCD connects the gauge-dependent gluon and quark degrees of 
freedom with the gauge-invariant properties of mesons. We study the p meson mass and decay 
constant for various versions of the gauge-dependent input functions discussed in the literature, 
which start to differ generically below the hadronic scale, and show qualitative different infrared 
behavior. We find that, once the gauge-dependent quark-gluon vertex is permitted to vary as well, 
the p mass and decay constant is reproduced equally well for all forms investigated. A possible 
conclusion from this is that these p-meson properties are only sensitive to changes in the input at 
scales above a few hundred MeV. 

PACS numbers: 12.38.Lg,14.40.Be,11.10.St 
Keywords: 



I. INTRODUCTION 

In the context of the Dyson-Schwinger equations 
(DSEs) of QCD, the homogeneous Bethe-Salpeter equa- 
tion (BSE) has been quite successfully used to study 
mesonic properties [l|-01, with a natural extension to 
baryons [5l-(7|. Such model calculations require both a 
numerical treatment as well as a truncation of the infinite 
coupled tower of DSEs, the latter of which may introduce 
artifacts in the results. Since the high-momentum behav- 
ior of correlation functions is essentially under control by 
perturbation theory, such artifacts will be most problem- 
atic at low momenta, i. e., at the order of the hadronic 
energy scale, or below. Unfortunately, this is also the en- 
ergy domain which is most important for hadronic prop- 
erties. Here, we investigate whether the whole low-energy 
domain is relevant, or only a part of it. The aim of this 
study is thus to investigate in particular to which extent 
the far infrared properties are of relevance. This inves- 
tigation is performed for the combined properties of the 
TT and p mesons, in particular, their masses and decay 
constants as a straight-forward and instructive indicator 
for the effects under consideration. 

The reason to suspect that only a certain energy win- 
dow, between a few hundred MeV and a few GeV, is 
relevant stems from two sources. On the one hand, in- 
dications have been found from various sources [8l4lOl| 
that at least for some mesons only the energy domain 
around a half to one fermi is significantly contributing 
to their properties. On the other hand, recent results 
on the gauge-dependent correlation functions of the ele- 
mentary constituents of mesons, the quarks and gluons, 
has given rise to the hypothesis that the far infrared be- 
havior may be qualitatively dominated by a gauge choice 
pT| - [l3| . as a consequence of the Gribov-Singer ambigu- 
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ity. This would again imply that the properties of mesons 
should be rather insensitive to the far infrared behavior, 
in particular of the elementary propagators and vertices. 

In fact, even the precise determination of such a simple 
correlation function as the gluon propagator for a single 
fixed gauge turned out to be a rather complicated prob- 
lem. Therefore, this whole problem set obtains also a 
different, purely technical, perspective: If the mesons are 
indeed rather insensitive to the far infrared, can then an 
incorrect input still yield the correct meson properties, 
and thus fool us into believing that we have obtained the 
correct input? 

For the present purpose two qualitatively different 
inputs are chosen, the so-called scaling-type [ll| and 
decoupling-type [ll|, [l4| - [l6j solutions for the Yang-Mills 
sector of Landau-gauge QCD. While the former shows 
an infrared power-law behavior, characterized by sc aling 
relations between the different correlation functions |17l - 
[131, the latter exhibit screening [ill, [l3-[l6|. At small 
volumes on a lattice, it is found that scaling-type and 
decoupling-type behavior can be associated with how 
the non-perturbative Gribov-Singer ambiguity is treated 
\\% . [l3l [21I [23 |. However, at the present time it is com- 
pletely unclear whether this persists in a continuum and 
infinite-volume formulation. 

The differences between both solution starts to man- 
ifest itself only at energies at or below a few hundred 
MeV, and becomes substantial only at even lower scales. 
Thus, they provide an ideal laboratory to study the ques- 
tions raised above. It is possible to take various different 
points of view with respect to these questions: 

1. If the lattice results [1, [0| are correct, it should 
not matter which of the solution is taken, since the 
p-meson's largest gauge-invariant contribution ap- 
pears to be only of the order of half a fermi large. 

2. Since it has been argued that only one of the so- 
lutions remains in the infinite- volume limit |14l4l6l | 
(which would then most likely be of decoupling- 
type [23,[2J|), it can be tested whether meson prop- 
erties would be at all able to distinguish between 



the solutions. In particular, since this implies at 
least one solution is wrong, this tells us whether 
meson properties are useful tools to identify errors 
made in the determination of the correlation func- 
tions in the deep infrared. 

3. Finally, the hypothesis that scaling and decoupling 
are just gauge choices can be tested with such a 
setup: If they are, the meson properties may not 
depend on the choice, though this requires to take 
into account the gauge-dependence of the quark- 
gluon vertex, like in QED [25|, [2a]. However, this 
cannot verify the hypothesis; it can at most falsify 
it. If the hypothesis would be assumed to be cor- 
rect, then the present study investigates how gauge- 
invariance is recovered in the Bethe-Salpeter equa- 
tion, and a lower bound for the gauge- dependence 
of the quark-gluon vertex is given. 

In any case, the main outcome of the present study will 
only be how sensitive some meson properties of some 
mesons are on the far infrared of the gauge-dependent 
correlation functions. 

For the present purpose, all model parameters are fixed 
to pion properties; then the p-meson is studied. Indeed, 
it turns out that the p is essentially not affected by the 
choice of Yang-Mills solution, as long as the quark-gluon 
vertex is adapted accordingly. This investigation com- 
plements similar ones, but with focus on different ob- 
servables [27r3l!l. 

The article is organized as follows: The input scaling- 
type and decoupling-type solutions are detailed in section 
inland the resulting quark propagator is studied in section 
mil The gauge-invariant results for the mesons are then 
presented in section IIVI A few concluding remarks can 
be found in section IVl 



II. RUNNING GAUGE COUPLING 

As noted above, two sets of solutions for the Yang- 
Mills sector will be used in the present study, both of 
them associated with the perturbative Landau gauge. 

One is the so-called decoupling-type solution [ll|, UA - 
lla . l23l . |24| . In this case, the gluon propagator is in- 
frared finite, thus exhibiting a screening mass, although 
it does not appear to have a pole mass [ll|, |32, ISj] • At 
the same time, the Faddeev-Popov ghost propagator is 
that of an essentially massless particle. Correspondingly, 
also higher correlation functions are expected to show 
a screening behavior although, except for some lattice 
studies in small volumes [3J|, this has not been analyzed 
in detail. 

The second solution is of the so-called scaling-type 
[nl \VA [l8j . It is characterized by correlation functions 
showing an infrared power- law behavior [13, H^ • ^1^ ^^~ 
ponents are tied to one single base exponent, namely the 
one of the ghost propagator [13, l20|- This propagator 




— sc 


- DCl 


- DC 2 


- ■ DC 3 



FIG. 1: The input running coupling ([T|. Here, and hereafter, 
SC denotes the scaling solution and DC the decoupling solu- 
tions. For more details, see 11]. 



has to be more divergent than that of a massless parti- 
cle. As a consequence, the gluon propagator has to be 
either infrared finite, like in the decoupling solution, or 
infrared vanishing, depending on the value of this critical 
exponent. 

As already mentioned, is has been proposed that both 
solutions are just non-perturbative gauge choices |lll - [l3| . 
but so far it has not be en p ossible to establish this beyond 
small lattice volumes |13l |. 

Irrespective of the solution, it is possible to define a 
running coupling by [111 ISSl j 



a{p')^a{^,'){G{p\^l')Yz{p\^l') 



(1) 



in the so-called miniMOM scheme [3a|. Herein G and 
Z are the dimensionless ghost and gluon dressing func- 
tions, respectively, which are obtained from the scalar 
part of the propagators by multiplication with p^. For 
the scaling-type solution this coupling is infrared finite, 
while for the decoupling-type solution it is infrared van- 
ishing like p^ . Note that this does not contradict the 
conjecture that both could be just gauge choices, since a 
definition of a running coupling in terms of ghost fields 
is by construction not gauge-invariant [33]. However, if 
preferred, it is possible to translate this running coupling 
into the standard MS scheme, at least perturbatively up 
to four loops [33. Nonetheless, for the present purpose, 
the form ([1]) is much more convenient. 

In the ladder truncation of the BSE coupled to the 
rainbow-truncated quark DSE, as we employ them here, 
this running coupling is actually the only way the Yang- 
Mills sector enters, and the individual gluon and ghost 
propagators are not needed. Here, the results for all solu- 
tions are taken from [ll|. Note that the decoupling- type 
solution is possibly not unique Jill, in agreement with 
small- volume lattice simulations [131]. For our purpose, 
we select three representatives, one of them (DC 3) be- 
ing in rather good quantitative agreement with lattice 
results [ll[ , in order to study also the dependence on dif- 
ferent variants of decoupling-type behavior. We fit these 
numerical results with a simple ansatz, which faithfully 
reproduces leading-order perturbation theory, along the 



lines of [37|. The resulting input functions are shown in 
figure [TJ Note that for simplicity we neglect the back- 
coupling of the quarks to the Yang-Mills sector, which is 
only a rather weak effect 33]. It can be seen that the 
coupling quantitatively depends on the choice of solution 
below a momentum scale of 1 GcV, and qualitatively be- 
low 100 MeV. 



III. QUARKS 

In this context, a meson is treated as a ouark-antiquark 
bound state via the homogeneous BSE [3l]. In such a 
setup, the dressed quark propagator S{p) appears as an 
input in the BSE. For a quark with momentum p, the 
propagator is given by 



z,(/) 



Sip) = Z,ip')- 






(2) 



with its two scalar dressing functions, the wave-function 
dressing Zq and the mass function M. To determine the 
quark propagator, we employ the rainbow truncation of 
the quark DSE. The quark DSE or QCD gap equation in 
its general form is given by 



Zmjn, 



s{p)-' 

S(p) = Zig^f D^,ip-q) 
J a 



Z2ii"f -p^ 



A" 



7p5(q)r^(p,q),(3) 



where J denotes the integral J d'^q/{2TiY regularized 
in a translationally invariant way with the regulariza- 
tion scale A, niq is the current-quark mass at the renor- 
malization scale, and the quark self energy E involves 
the renormalized dressed gluon propagator D^^ and the 
renormalized dressed quark-gluon vertex r°(p, q). Here, 
a is a color index. A" the Gell-Mann SU{2>) color matri- 
ces, g the strong coupling constant, and the Zi are renor- 
malization constants (for more details and a full account 
of the renormalization procedure, see e.g. |39|). 

In rainbow truncation in Landau gauge the contribu- 
tions from the dressed gluon propagator and the dressed 
quark-gluon vertex are contracted to a single, effective 
object, for which we take the ansatz 
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a(fc2) / „ „£) 









(4) 
(5) 



The particular functional form of the effective interaction 
chosen here is inspired by various model ansaetze used in 
the past which for the most part make use of the knowl- 
edge of the perturbative QCD running coupling as well as 
a simple parametrizations of the low- and intermediate- 
momentum parts of the effective interaction |39l - l43| . Also 
note that this has been done before in a similar way in 
\4M. albeit with a different ansatz for the vertex. 
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FIG. 2: The quark wave-function dressing (top) and mass 
function (bottom) for the four choices of parameters of the 
quark-gluon vertex fixed to pion properties, see table [III 



The parameters D and uj are fixed by the pion prop- 
erties as discussed below. At the same time, one has to 
fix a value for the current quark mass which acts as a 
driving term for the non-perturbative M in the quark 
DSE. Specified in analogy to Ref. [4^ we have a value 
of rtiq = 0.0047 GeV at our renormalization point of 
(19 GeV)2 and we use two degenerate fiavors. The result- 
ing quark dressing functions are shown in figure [5J They 
no longer show any qualitative difference when compar- 
ing results corresponding to the scaling- and decoupling- 
type inputs; in fact, only a small quantitative difference 
remains. The qualitative uniformity is a direct conse- 
quence of dynamical chiral symmetry breaking, since it 
makes the quark effectively decouple in the infrared and 
thus blind to the qualitative differences. The smallness 
of the quantitative differences among the solutions for 
Zq and M turns out to stem from the requirement of 
accurately reproducing the pion properties. 

Note that some results are already available for the 
quark-gluon vertex and its further tensor structures as 
well as its impact on bound-state calculations (see e.g. 
[45l44g | and references therein). However, these are either 
specific to the scaling case [45| , from the lattice for small 
volumes [46l |. or from more involved calculations using 
functional methods [47h49| and cannot be used in the 
same way as the gluon input for the comparison aimed 



at here. 



IV. MESONS 

In our study of meson properties, we employ the ho- 
mogeneous quark-antiquark BSE, which in general form 
reads 



r(fc;P) 



Kik;q-P)Siq+)T{q;P)S{q-), (6) 





m^ 


u 


nip 


fp 


Scaling 


0.1417 


0.0932 


0.7443 


0.1492 


Decoupling 1 
Decoupling 2 
Decoupling 3 


0.1420 
0.1401 
0.1389 


0.0956 
0.0897 
0.0805 


0.7641 
0.7155 
0.6321 


0.1533 
0.1439 
0.1239 


Experiment 


0.1396 


0.0922 


0.7755 


0.1527 



TABLE I: Results for tt and p mass and decay constant after 
adjusting the parameters of the vertex ansatz (|4]) such that 
the scaling solution best reproduces the pion properties. All 
quantities are given in units of GeV; our results are accurate 
to less than half an MeV. 



where q+ and q_ are the quark- and antiquark momenta, 
respectively. The Bethe-Salpeter amplitude r(q; P) de- 
pends on their total and relative momenta P and q with 
the relations q^ = q + rjP and q- — q — {1 — ri)P. The 
homogeneous BSE is only valid on-shell, i.e. the total 
momentum of the bound state satisfies the (Euclidean) 
on-shell condition P^ — —IVP, where M is the mass of 
the meson under consideration, r/ with < 77 < 1 is the 
momentum partitioning parameter, which is arbitrary in 
the Lorentz-covariant context. The 2-fermion irreducible 
quark-antiquark interaction kernel K in general depends 
on the relative momenta as well as the total momentum 
of the meson. In ladder truncation, the kernel does not 
depend on the total momentum, but only on the differ- 
ence in relative momentum (or in other words on the 
gluon momentum) k — q, and one obtains 



r(fc;P) 
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(7) 



Note that the Dirac structure hidden in Eq. ^ has been 
made explicit through the Dirac 7 matrices and the or- 
dering of the terms under the integral has been chosen 
such that Dirac indices can be omitted. 7"^ denotes the 
transversely projected j^-jriK - qT){k^ - q^)/{k - qY, 
and a and F are the same as in Eq. ([S]). 

As becomes clear from Eq. ([6]), the input quark prop- 
agator is needed for the arguments g±, which are com- 
plex vectors due to the on-shell condition. Therefore, 
an analytic continuation of the propagator is necessary. 
In this context, a complication arises in the case of the 
scaling-type coupling. Since it does not go to zero as 
p^ — >■ 0, it leads to a singularity in the vector part -g- 

of the inverse quark propagator [50, [5l|. Thus, in order 
to treat scaling- and decoupling-type inputs on the same 
footing, we employ the method detailed in f52i| to obtain 
the solution of the gap equation in the complex plane. 
The method has two main features: First, instead of the 
quark momentum as in Eq. ([3]), the gluon momentum 
is used as integration variable. As a consequence, the 
coupling is only needed for real and positive arguments, 
whereas during the iterative solution the quark propaga- 
tor itself has to be evaluated for complex momenta, if the 
external momentum p is complex. This, however, leads 
to the second feature: By demanding that J- and 



M_ 



are analytic in the domain of interest, the gap equation 
can be iterated on the boundary of that region, and the 
values inside this boundary can be obtained in a numer- 
ically stable way from an appropriate implementation of 
the Cauchy formula [5j] . During the iteration, all singu- 
larities are therefore eliminated and the result represents 
a solution of Eq. ([3]) in the complex plane which by con- 
struction does not exhibit the kinematic singularity in 
the case of the scaling-type input, such that a consistent 
treatment of both types of input is possible. Once the 
quark propagator is provided in such a way, the homo- 
geneous BSE can be solved numerically (for more details 
on the techniques used in this context see j54|). 

Together with the rainbow truncation of the quark 
DSE, the ladder truncation of the BSE is used, since this 
so-called rainbow-ladder truncation satisfies the axial- 
vector Ward-Takahashi identity in QCD, which encodes 
chiral symmetry and related properties of the theory in 
the context of the DSEs |55| . Satisfaction of this identity 
in a model setup such as the one used here guarantees the 
correct implementation of chiral symmetry and its dy- 
namical breaking and produces, e.g., the well-known re- 
sult of a massless pion in the chiral limit. Via the identity 
one finds a generalized version of the Gell-Mann-Oakes- 
Renner relation that is valid for every pseudoscalar me- 
son, see [56| and references therein. With the pion an- 
chored in this way in such a calculation, its (and also 
kaon) properties have been subject to successful study 
for the past years, exemplified by the excellent descrip- 
tion of the pion's electromagnetic properties, e.g. |57l.l58|. 
With similar success, vector meson properties have been 
investigated, e.g. decay constants, hadronic decay widths, 
as well as electromagnetic form factors and charge radii 

Thus, the tt and p are ideal objects for the purpose of 
the present study. We present results for masses and 
decay constants, where the latter are sensitive to the 
Bethe-Salpeter amplitude of the meson under considera- 
tion. The numbers are compiled in two tables to highlight 
two different aspects of the investigation. First, we keep 
the same vertex ansatz for all four different Yang- Mills in- 
put sets. We adjust the parameters of the vertex ansatz 
dH) such that the scaling solution reproduces the pion 
mass and decay constant adequately. The resulting val- 
ues for D and w are 0.8 GeV* and 0.8 GeV, respectively. 
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for the optimized fits. More precisely, we investigate the 
combined effect of the couphng and vertex dressing func- 
tions on different momentum scales. This is shown in 
figure El It is nicely visible that the four cases for the 
ansatz start to differ appreciably only below 1 GeV, but 
are qualitatively different in the far infrared. This ob- 
servation supports the argument that the p has a gauge- 
invariant structure which exists at a level of 0.5 fm or 
below [8, 9], and is not too sensitive to larger scales. 
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FIG. 3: The momentum dependence of the ansatz Q for 
the parameter values D and uj of table |I] (top) and for the 
optimized values of table |TI] (bottom). 



Then we compute the masses m and decay constants /, 
"^tt; fir, 'iTT'p, and fp with the same parameters also for 
the different decoupling-type cases; the corresponding re- 
sults are listed in table U in units of GeV. One can see 
that the more the input coupling differs from the scaling 
case, the larger the differences for /tt, nip, and fp become, 
while the general picture is not spoiled at a qualitative 
level. The small variation in ttt,^ is due to the infiuence of 
the axial- vector Ward-Takahashi identity. This increase 
in discrepancy when only changing the gluonic input has 
also been observed in the literature [30, |311 • 

Next, we offer another possibility: irrespective of 
whether the different inputs correspond to different 
gauges or are just different truncation assumptions, it is 
to be expected that the quark-gluon vertex should change 
accordingly when the running coupling changes. To in- 
vestigate this scenario, we have fitted the parameters D 
and uj for each input individually to reproduce the pion 
parameters. The corresponding results are shown in ta- 
ble ini As one can clearly see, already a very moderate 
change in the parameters without excessive finetuning 
yields agreement of rrijr, J-k, nip, and fp for all inputs 
on the percent level. Thus, any consistent change of 
the quark-gluon vertex and the gluonic sector permits an 



equally good description of m.^-, /^, 



i-p. 



and fp. While 



not guaranteed, it is at least possible that a similar state- 
ment holds for arbitrary other observables as well; some 
confirmation of our observation along these lines has al- 
ready been found in [23, [2^ . 

It is interesting to check how the full ansatz (U) looks 



V. CONCLUSIONS 

Summarizing, we have shown that in the coupled sys- 
tem of quark DSE and meson BSE in rainbow-ladder 
truncation, changing the gluonic input in these equa- 
tions can be compensated, at the level of the gauge- 
invariant hadronic properties, by changing the corre- 
sponding quark-gluon vertex accordingly. This becomes 
apparent only at the level of the gauge-invariant bound- 
state properties, in our case meson masses and decay 
constants, but not on the level of the gauge-dependent 
quark propagator. 

Returning to the three different questions posed in the 
introduction, this pattern in our results can be inter- 
preted as follows. For this two central observations are 
important. First, while the meson properties are qualita- 
tively unchanged when changing some part of the input, 
there is some sensitivity to the far infrared left. Other- 
wise the results in tableUwould agree. Thus, the infrared 
is not irrelevant, and its behavior cannot be chosen arbi- 
trarily. However, if all inputs are varied consistently, the 
meson properties are unchanged, see table |lTl despite the 
fact that the infrared is still qualitatively different for the 
inputs, as is visible in figure[21 The reason for this insen- 
sitivity is found in the quark propagator, as discussed in 
section Hill Chiral symmetry breaking screens the quark 
from the low-energy interactions, and thus makes it blind 
to the far infrared. The conclusions may thus be altered 
in the absence of chiral symmetry breaking, but then, the 
set of possible solutions in the gluon sector may differ as 
well. An indication of such a change is found in models 
with scalar instead of fermionic matter [60| . 

This leads to the following conclusions: 

1. From a technical point of view, it is not possible to 
deduce the infrared behavior of the input gauge- 
dependent correlation functions from the meson 
properties alone. Reproducing meson properties 
therefore not necessarily implies a faithful solution 
for the gauge-dependent correlation functions. 

2. Though rather insensitive to the infrared, the p me- 
son is still infiuenced by the long-distance physics. 
Only a consistent long-distance behavior can repro- 
duce the meson spectrum adequately. 

3. If scaling and decoupling are indeed mere gauge 
choices, then, not surprisingly given the experience 





D 


ijj 


m^ 


h 


nip 


/p 


Scaling 


0.800 


0.800 


0.1417 


0.0932 


0.7443 


0.1492 


Decoupling 1 
Decoupling 2 
Decoupling 3 


0.790 
0.820 
0.870 


0.810 
0.785 
0.760 


0.1414 
0.1411 
0.1407 


0.0932 
0.0936 
0.0935 


0.7443 
0.7484 
0.7485 


0.1496 
0.1502 
0.1514 



TABLE II: As tableU] but D and cj have now been adjusted for each input separately to fit the pion's mass and decay constant. 
All quantities are given in units of GeV except for D which has units of GeV^. 



of QED [25|, |26| , there is a gauge-dependence also 
in the matter sector. If this is taken into account, 
the results for mesons are gauge-invariant. How- 
ever, in agreement with direct investigations |60| . 
the gauge-dependence in the matter sector is rather 
small. 

Thus in summary, the infrared region requires a con- 
sistent treatment when quantitative precise results are 
aimed at, but there is more than one infrared behavior 
leading to the same results. This is also an important 
result for the eventual aim of a self-consistent inclusion 
of the quark-gluon vertex in the future: Truncation ar- 
tifacts in the far infrared are of smaller relevance than 
those at mid-momentum. 

Regardless of this interpretation, we have reported 
results which can be taken as a guiding line for how to 
proceed when aiming at getting gauge-invariant results 
from gauge-dependent correlation functions, as sup- 
ported also by other investigations [28]. For this study. 



we have chosen a straight-forward, but limited set of 
gauge-invariant observables, and we have made this clear 
together with the limits of the interpretations provided, 
as appropriate. Naturally, further investigations along 
these or similar lines must follow to shed more light on 
this important set of problems, which could be along 
the lines of corresponding QED investigations [25|, |26| . 
A particular challenging and interesting topic would be 
Regge trajectories. 
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